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ABSTRACT 


The  two  dimenelonal  stendy-stftte  probleo  of  the  effect  of  a  step  pressure 
traveling  with  superseismic  velocity  on  the  surface  of  a  half-space  is  treated 
for  an  elastic -plastic  material.  The  plasticity  condition  selected  is  suitable 
for  a  granular  niedlim  where  inelastic  deformations  are  due  to  internal  slip 
subject  to  Cbulomb  friction. 

The  problem  is  inherently  nonlinear  and  leads  to  a  system  of  coupled 
differential  equations  which  are  solved  by  digital  cornputer.  Numerical 
solutions  are  tabulated  as  functions  of  the  significant  nondlmenslonal  parameters, 
i.e.  of  the  Mach  nuitbeir,  Poisson’s  ratio  and  of  a  value  a  defining  the  internal 


friction. 
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I  INTRODUCTION 


The  two  dimensional  problem  of  the  effect  of  a  pressure  pulse  p(x  -  Vt) 
progressing  with  the  velocity  V  on  the  surface  of  an  elastic  half-space. 

Fig.  1,  has  been  treated  by  Cole  and  Huth  [1]  for  a  line  load  and,  by  super¬ 
position,  may  be  found  for  any  other  distribution  p(x  -  Vt).  Miles  [2]  has 
considered  the  three  dimensional  problem  of  loads  with  axially  symmetric 
distribution  p(r,t)  over  an  expanding  circular  area  on  the  surface,  Fig.  2. 

He  has  demonstrated  that  the  plana  problem  [1]  contains  the  as^tiptotic  so¬ 
lution  for  the  three  dimensional  problem  [2]  in  the  region  near  the  wave  front. 
The  actual  solution  of  the  three  dimensional  problem  would  require  a  great 
numerical  effort, which  can  be  avoided  by  using  the  solution  of  the  plane 
problem  to  estimate  the  effect  of  circularly  ejqpanding  surface  loads. 

Real  materials  can  not  be  expected  to  be  elastic,  and  solutions  of  the 
three  dimensional  problem.  Fig.  2,  for  dissipative  materials  are  hopelessly 
complex.  However,  estimates  for  the  three  dimensional  case  can  be  made  from 
generalizations  of  the  problem  treated  in  [1]  for  dissipative  materials. 

This  has  been  done  for  linearly  viscoelastic  materials  by  Sackman  [3],  and 
Workman  and  Bleich  [U],  in  the  superseismlc  and  subseisraic  ranges,  respectively. 
For  possible  application  to  granular  media  the  present  report  considers  an 
alternative  material  where  internal  slip  subject  to  Coulomb  friction  may  occur. 
The  problem  has  previously  been  considered  by  Bleich  and  Heer  [5J  for  the 
range  of  low  subseismic  velocities  V,  while  the  more  interesting  superseismic 
case  is  the  concern  of  this  report. 

The  slip  mechanism  in  the  medium  makes  the  problem  nonlinear,  such  that 
superposition  is  not  permitted  and  each  pressure  distribution  p(x  -  Vt)  poses 
a  separate  problem.  The  present  report  treats  the  case  of  a  progressing  step 
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load  p(x  -  Vt)  ■  p^H(Vt  -  x).  An  approach  permitting  a  simplified  solution  of 
the  important,  but  very  complex  case  of  a  decaying  surface  pressure  is  discussed 
in  Section  IV. 


Based  on  concepts  of  the  theory  of  elastic-plastic  materials,  Drucker  and 
Prager  [6]  have  shown  that  a  material  subject  to  internal  Coulomb  friction  can 
be  represented  by  an  ideal  material,  the  behavior  of  which  is  governed  by  a 


plastic  potential 


F  -  l/Jgl  + 


where  and  are  the  invariants 


h  ‘  "li 


while  or  >  0  and  k  >  0  are  material  constants,  a  is  related  to  the  angle  of 
internal  friction  and  is  therefore  subject  to  the  limit  a  i  ,  [6],  and  k  is 
a  measure  of  the  cohesion.  Because  the  surface  pressures  for  which  this  study 
is  intended  are  large  compared  to  the  magnitude  of  cohesion,  it  suffices  to 
consider  the  limit  k  0,  giving  the  simpler  plastic  potential 

F  =  l/Jgi  +  crJ^  .  (4) 

The  behavior  of  the  material  is  described  by  the  following  statements: 


1.  To  represent  a  granular  material  with  no,  or  very  small  cohesion,  the  mean 
stress  -j-  must  be  compressive,  or 

Jj^<0  .  (5) 


2.  If,  in  an  element  of  the  raaterial  at  a  given  instant. 


4 


P  <  0 


(6) 


the  changea  in  atreas  and  atrain,©^^  ,  Vijl  ha  related  hy  the 
conventional  elastic  relations « 


3>  Hoeever^  if  the  yield  condition  at  a  tiiae  t  ia  satisfied 

T  m  0  (7) 

three  possibilities  exist.  There  siay  be  fiurtber  loading  of  the  element 
with  permanent  plastic  deformation  and  dissipation  of  energy^  in  which 
case  F  ■  0.  Alternatively,  there  may  be  unloading  without  permanent 
deformation,  in  which  case  F  <  0.  Finally,  there  »ay  be  a  neutral  state 

e 

where  F  ■  0,  but  without  permanent  deformation  or  energy  dissipation. 

For  the  first  case,  with  plastic  deformation^  the  total  strain  rate  will 
be  the  sum  of  an  elastic  and  a  plastic  portion 


•  •£  •? 
•ij  -  ‘ij  ^  ‘ij 


(8) 


where  is  obtained  from  the  conventional  elastic  relations,  while 


(9) 


X,  ^ich  must  be  positive. 


X  >  0 


(10) 


is  an  a  priori  unknown  function  of  space  and  tine. 
In  case  of  unloading,  and  in  the  neutral  case  the 
relAtioDi 


« 


•ij 


(11) 


app]^>  The  neutral  case  occurs  in  the  solutions  Obtained  in  regions 
without  change  in  stress  or  strain. 


5 


The  fact  that  the  same  set  of  differential  equations  does  not  hold  everywhere, 
but  that  there  are  regions  with  moving,  a  priori  unknown  boundaries,  complicates 
the  solution  of  dynamic  problems  in  this  type  of  material  considerably.  In  the 
following,  the  basic  equations  will  be  formulated  separately  in  regions  with  and 
without  permanent  deformations  at  the  particular  time  t,  and  the  solutions  will 
be  matched  to  obtain  a  complete  solution  satisfying  the  prescribed  surface  con™ 
ditionso  The  problem  being  much  too  complex  to  expect  closed  solutions,  a  numeri¬ 
cal  approach  suitable  for  digital  computers  will  be  employed.  The  technique  is 
related  to  the  theory  of  characteristics  and  is  a  generalization  of  the  method 
used  by  Bleich  and  Nelson  [?]« 

The  problem  to  be  solved  considers  only  the  steady-state,  i.e.  the  fact  is 

ignored  that  in  reality  the  loads  p(x  -  Vt)  in  Fig.  1  must  have  begun  at  some 

large  but  finite  negative  value  of  time.  This  omission  of  the  initial  condition 

results  in  a  lack  of  uniqueness,  which  can  be  removed  by  consideration  of  the 

character  of  solutions  of  the  problem  in  Fig.  2.  The  lack  of  uniqueness  and  the 

remedial  consideration  is  best  seen  in  the  elementary  exaji5)le  of  a  half-space  of 

an  inviscid  con^resslble  fluid  loaded  by  a  uniform  pressure  pulse,  p,  which 

progresses  with  supersonic  velocity,  V  >  c.  There  is  an  obvious  solution, 

Fig.  3a-,  in  which  the  load  produces  a  plane  wave  of  intensity  p  progressing  wltn 

-1  c 

a  front  inclined  at  the  appropriate  angle  f  «  sin  -  .  However,  this  is  not  the 
only  steady-state  solution.  An  alternative  is  a  plane  wave,  the  front  of  which  is 
inclined  at  the  angle  l8o°  -  f  .  Combinations  of  the  two  solutions  are  also 
correct  steady-state  solutions.  To  find  states  generated  by  the  application  of 
pressure  on  the  sturface  only,  it  can  be  reasoned  that  solutions  which  include  the 
wave  front  shown  in  Fig.  3t>  can  not  apply  because  the  medium  ahead  of  the  front 
shown  in  Fig.  3a  should  be  vuidisturbed  when  the  applied  load  moves  with  super¬ 
sonic  /elocity.  Thus,  in  case  of  the  fluid  a  unique  solution  is  obtained. 

Similar  reasoning  will  be  used  in  the  body  of  the  report. 
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II  FORMULATION  OF  TIE  BASIC  B^UATIONS 


Figure  4  indicates  the  half»flpace  and  a  system  of  Cartesian  coordinates. 

X  is  in  the  direction  of  motion  of  the  step  load,  y  and  z  are  normal  to  the 
surface  in  and  out  of  the  plane  of  the  paper,  respectively.  The  analysis  con¬ 
siders  the  case  of  plane  strain,  t  ■  0,  when  the  velocity  V  of  the  step  load  is 
superseismic,  i.e.  lai'ger  than  the  velocity  of  P  waves  in  the  material  when  slip 
does  not  occur.  Throughout  the  analysis  it  is  assumed  that  the  strains  are  small. 

As  stated  in  the  introduction  there  are,  in  general,  inelastic  regions  in 
space-time  where  permanent  deformations  with  energy  dissipation  occur,  and  other 
regions  where  changes  of  stress  and  strain  are  entirely  elastic.  The  basic 
differential  equations  for  the  two  types  of  region  must  be  treated  separately. 

In  addition  it  will  be  necessary  to  consider  the  possibility  of  shock  fronts, 
i.e.  degenerate  infinitely  narrow  regions  where  the  differential  equations  break 
down. 


a.  Inelastic  regions 


Combining  the  familiar  elastic  stress  strain  relation 

1+v 


.E 

•ij 


“ 

0 


^ij  “  1+v  ®ij  ^kkj 


(12) 


where  6, ,  is  the  Ktoneker  delta,  and  the  relation 


‘iJ 


(13) 


expressing  the  strain  rates  in  terms  of  the  velocity  components,  gives,  for 
plane  strain,  four  constitutive  equations 


-v(d  +6). 
ax  E  L  ^  '  yy  zz'j 


+  \ 


JL 

ao 


(14) 


XX 
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F  I  6  5b 
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0  »  - 
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,  a  -  V  (0  +  0  ) 

L  JO'  XX  zz  , 


^  -  V  (0  +  ) 

zz  '  yy  XX  . 


^  1  i  ^  ^  aF 

~  +  ~a—  f  +  X  — 

?»y  ^x  G  hr 


30 

(15) 

yy 

30 

(16) 

zz 

(17) 

where  u  and  v  are,  respectively,  the  x  and  y  components  of  the  particle  velocity. 
Further,  there  are  two  equations  of  motion,  which  are  in  a  linear  theory 


ax  ay  ”  ^  at 


(18) 


|L  +  . .  5i 

ax  ay  ^  at 


(19) 


The  yield  condition,  Eq.  (U),  and  Eqs.  (14)- (19)  form  a  set  of  seven 
equations  governing  inelastic  regions.  Inherently,  however,  the  unknown 
function  X  must  satisfy  the  inequality  X  >  0,  Eq.  (10),  required  for  an 
element  in  an  inelastic  region. 


It  was  found  convenient  to  express  the  four  unknown  stresses  by  four 
other  variables:  the  invariant  ,  the  two  principal  stress  deviators  s^ 
and  Sg  and  the  angle  9  formed  by  the  direction  of  s^^  with  the  surface,  Fig.  4. 
The  appropriate  relations  are 


2  2  1 

^xx  "  ®2  6  +  cos  0  +  j 

2  2  1 

0yy  »  8^  sin  8  +  Sg  cos  9  +  j 


0  »  -  8,  -  +  -  J, 

zz  1231 


T  ■ 


01  -  02 


sin  20 


(20) 

(21) 

(22) 

(?3) 
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In  the  nvunerical  analysis  the  subscripts  1  and  2  vill  be  selected  such 


that  is  the  roa,;or  compressive  deviator. 


Because  the  steady  state  case  is  considered,  all  quantities  appearing  in 
the  analysis  which  are  functions  of  x  and  t  must  be  of  the  form  f(x  -  Vt).  For 
the  step  load  p  =  P^HCVt  -  x),  dimensional  considerations  similar  to  those  used 
in  [5j  make  it  plausible  that  the  various  quantities  do  not  depend  on  x  -  Vt  and 
y  separately,  but  must  be  solely  functions  of  the  variable 


,  X  -  Vt 

§  M  - 

y 

or,  alternatively,  of  the  angle  qp,  shown  in  Fig.  4,  and  defined  by 


(24) 


5  ■  cot  . 


(25) 


The  transformation  to  the  new  independent  variable  f ,  which  will  be  seen 

to  be  successful  in  obtaining  a  solution,  changes  the  partial  differential 

equations  obtained  above  into  a  set  of  simultaneous  ordineiry  differential 

equations.  Noting  ^  ^ —  , 

^  sin*^  f 


p 

-?L  -  i  A  a  »  sin'  tp  d 

5x  ”  y  d5  “  "  '  y  d(p 

i  5  d  sin  2o  d 

ay  y  d5  2y  dip 

B  V  d  V  .2  ^  d 

at*’‘y^"y® 


(26) 


and  defining 


L-+ - (27) 

crJ^  sin  «p 

the  seven  basic  equations  become,  respectively, 


*1  ®1®2  ^2  " 


(28) 
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roiM 


0 


(29) 


*  >11 
cos  0  sin  (<p  -  0)  +  sin  6  cos  (tp  -  0)  s^  +  sin  ®P  ^  " 

-  cos  (qp  -  20)(s^  -  Sg)  0  +  pV  sin  cp  u  =« 

i  *  1  * 

sin  0  sin  ((p  -  9)  s^  -  cos  0  cos  (cp  -  9)  Sg  -  cos  cp  ^  + 

+  sin  (<p  -  20)(Sj^  -  Sg)  0  +  pV  sin  9  V  =  0 


2  •  2  '  l-2\)  1  ’  X  * 

cos  0  s,  +  sin  0  8_  +  r -  -r  J,  -  sin  26  (s,  -  s-)  6 

1  2  1+v  31  1  2 

-  [s^  cos^  8  +  Sg  sin^  0  -  ^  ^  f  °  ^31) 


2  '  2  •  l-2\)  1  •  V  ' 

sin  0  s-  +  cos  0  s-  +  r — =“  ■=•  J,  +  sin  29  (s^  -  s.)  0 

1  2  1+v  31  1  2 

-  [s^  sin^  ®  ®2  ®  ~  ~  cot  9  v'  =  0  (32 

>11  >1  nr 

sin  20  s^  -  ^  sin  20  Sg  +  cos  20  (s^  -  Sg)  0  “  f  (®i  "  Sg)  sin  20  ^  - 

-|cot9u  +|v  =  0  (33) 


l"2v  Cx.2t  GL  -  G  2G  ^ 

ITT  *^1  -  T  "  ^  V  +  Y*  cot  9  V  =0 


where  primes  indicate  the  derivative  with  respect  to  9. 

Differentiation  of  Eq.  (28)  yields  a  seventh  differential  equation 

(2si  +  Sg)  s^  +  (2Sg  +  s^)  Sg  -  =  0  (3!?) 

Eliminating  the  velocities  u  and  v  from  Eqs.  (29)- (3^)  leaves  a  set  of 
five  differential  equations  in  five  unknowns  which  are  related  to  the  stress 


pattern  alone.  Defining,  for  convenience,  the  angle  y  (see  Fig.  k)  and  the 
quantity  X  ■  X(ep) 


(32) 


Y  -  -  e 


(36) 


X  - 


.  2 

sin  9 


(37) 


the  five  equations  become: 
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-  6a^J^ 

0 

0 

GL 

V 

_ 

(38) 


This  set  of  equations  is  linear  and  homogeneous  in  the  four  derivatives  and  in 
the  quantity  L,  and  may  be  satisfied  by 

61  =  Sg  =  Ji  »  (Si  -  Sg)  0  »  L  «  0  (39) 

However,  L  »  0  implies  X  «  0,  which  violates  Eq.  (10).  It  follows  that  in 
an  inelastic  region  the  determinant  of  Eqi.(38)  must  vanish,  giving  the 
"deterrainantal  equation" 


^  V3 


0 


(40) 


where 
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and 


«  2  [1  +  (1  -  2X)(1  -  2v)J 
bg  »  9  C08  2y  +  (1  -  2X)  [1  -  2v  “  (l+v)i 

b^  »  (1  +  v)(l  -  2X)(1  +  2\^f  -  a^X 


- 


or^J, 


®2 


/ 


(41) 


(42) 

(43) 


Due  to  the  vanishing  of  its  determinant  only  four  of  the  five  Eqs»  (38) 

tit  I 

are  independent.  As  L  may  not  vanish,  s^  ,  s^  ,  and  9  can  always  be 
expressed  in  terms  of  L, 


^1*1  *’4^  T 

'  1  /  .  ^/v  V  ^  GL 

=2  "  2  (=1  *  "2><^  -  T 

s'  =  *>6  f 

Jl  =  3(s^  +  Sg)  b^  ~ 

Velocities  and  accelerations,  if  desired,  are 


(44) 

(45) 

(46) 
(4?) 


u  = 


2pV 


OL 


[^4  sin(2Y  -  9)  -  Sgbg  cos(2y  -  cp)  »-  (b^+  2b^)  sin  9)  j 


(48) 


I  “  (S-l"^  ^  (ri 

^pvTrn"^  ^^^4  -  9)  +  29bg  sin(2Y  -  9)  -  (^5+  2b.^)  cos  9]  ^ 


V  = 


(4‘^) 


••  V  ,  2 

u  =  —  sin  9  u 

y 

••  V  ,  2  .• 

V  =  —  sin  9  V 

y 


(50) 

(51) 
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where 


4  “  1-2X 


cos  2y  -  3X 


»  I  (1+v) 


b 

l-2v  ^ 
1+v  b- 


«  y ; 


+  l+  crV3  +  p 


,  sin  2y  2 

'^e'fTrlxT  ^ 
”7  ■  ^  ^ 


(52) 


Since  Eq.  (40)  must  remain  valid  throughout  an  inelastic  region,  it  may 
be  differentiated  with  respect  to  cp*  This  leads  to  an  expression  which  contains 
the  first  derivatives  of  the  stresses  linearly,  so  that  substitution  of  Eqs. 
(44)-(47)  into  this  ejqpression  furnishes  a  linear  equation  for  the  value  of  L; 


'  p 

2Xb  (l-2v)  sin  2f  +  ^  Xb  sin  2^  [2(l+v)(l+^) 

+  3^]  +  ] 

+  2b«[3  sin  2y  sin^  9  +  X  sin  2np  (l-2v)  -•  4X  sin  29  (1+v) 

k  ^ 

b^[(l+v)(l+2U)(l-2X)(b^+  -  Xbi^p]  + 

<  +  bgCbi^  cos  2y  +  20bg  sin  2y  +  b^(l-2X)(l-2v)  - 

,  2 

'  sin  cp 

-  12b  or^(l+v)(l-2X)] 

> 

(53) 


III  I 

The  derivatives  3^,82,  and  6  can  be  obtained  by  substitution  of  Eq.  (53) 
into  Eqs.  (44)-(47). 


If  the  values  of  8^,8^,  and  6  are  known  on  one  boundary  of  an 
inelastic  region,  their  values  in  the  interior  of  such  a  region  can  be  found 
by  forward  integration.  Of  course,  the  starting  values  of  s^  ,  Sg  ,  and  0 
must  satisfy  the  yield  condition,  Eq.  (28),  and  the  determinantal  equation, 
Eq.  (4o).  Further,  the  integration  can  be  carried  on  only  when 
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L  <  0 


which  condition  follows  from  Sq.  (2?)  noting  X  >  0,  Eq.  (10),  and  <  0. 


Inelastic  shock  fronts 


The  analysis  of  the  previous  subsection  treated  regions  of  finite  extent, 
and  the  additional  possibility  of  infinitely  thin  regions,  i.e.  shock  fronts, 
remains  to  be  considered.  If  such  fronts  exist  in  the  present  problem  the 
equations  obtained  above  should  indicate  this  by  becoming  singular,  since  at 
least  one  of  the  derivatives  of  the  stresses  must  become  infinite  at  a  shock 
front.  Instead  of  searching  for  singularities  it  is  better  for  a  physical 
understanding  to  demonstrate  the  existence  and  properties  of  shock  fronts  in 
general.  This  general  derivation  automatically  answers  the  question  of  stability 
of  the  fronts,  by  proving  that  a  front,  where  the  stress  rises  with  an  arbitrarily 
steep  slope.  Fig.  5a;  will  not  disperse  but  propagate  without  change  of  slope. 

Consider  the  basic  Eqs.  (l4)-(19),  which  apply  to  any  type  of  wave 
propagation  in  plane  strain.  To  investigate  the  possibility  of  plane  pressure 
waves  without  shear,  the  y  direction  is  selected  as  the  direction  of  propagation. 
For  such  a  wave,  the  shear  t,  the  hoi'izontal  velocity  u  and  all  derivatives  with 
respect  to  x  vanish,  while  for  reasons  of  symmetry  a  ■  o  .  Since  x,  y,  z  are 

the  principal  directions,  the  stress  deviators  s  and  s  become  s  =  s  ,  s  =  . 

X  y  X  c  y  1. 

using  the  relation  Sg  *  -  J  s^  ,  the  yield  condition  (28)  becomes 


while  Eqs.  (l4)-(19)  furnish  three  independent  relations 


ay  2G  ®1  9K  '^1  ^  ^®1  "  ^  “^l^ 
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(57) 


where 


1  . 

5g 


1.  j 

9K 


2a^J^] 


(58) 


(59) 

(60) 


If  a  front,  at  which  inelastic  deformation  occurs,  is  to  propagate  with  a 
velocity  c  without  change  of  slope,  it  is  necessai’y  that 


a^f  (y  -  ct)  =  a^f  (C ) 

(61) 

agf  (y  -  ct)  =  a^f  (C) 

(62) 

a^f  (y  -  ct)  -  a^f  (C) 

(63) 

where  C  “  Y  ”  ct,  the  values  of  a^  and  c  are  free  constants,  while  f  is  an 
arbitrary  function.  Since  the  signs  of  the  coefficients  are  undefined,  one 
can  select  the  sign  of  f,  considering  fronts  of  the  type  shown  in  Fig.  5a. 
Selecting 


f  >  0  (64) 

4*V»rt  ^  I  A 

v..cirOv.  wx  1^1 

f’  <  0  (65) 

where  the  symbol  '  indicates  derivatives  with  respect  to  Q, 


IT 


j 

Substituting  Eqs.  (6l)-(63)  into  Eqs.  (55)- (58)  and  eliminating  the 
velocity  v  yields 


0 


(66) 


f 


1 

2 


1  -2 
i  fi£- 
3  3K 


5 


pc  Xf 


I  aj^  + 


=  0 


(6?) 


f 


2 


^2  2 
0'a2 


(68) 


If  shock  fronts  of  the  type  sought  exist,  these  three  equations  must  be 
satisfied  for  arbitrary  fvinctions  f,  subject  to  the  limitations  of  Eqs.  (61|), 
(65)  and  subject  to  the  condition  X  >  0.  Equation  (68)  gives 


(69) 


Equations  (66)  and  (67)  permit  nonvanishing  values  f,  X  and  f*  only  if  the 
determinant  of  the  coefficients  of  f ’  and  pc  Xf  vanishes,  yielding  after 
substit’ition  of  Eq,  (69) 


-2 

c 


K 

P 


(1  +  20/3  f 


(70) 


However,  the  result  is  valid  if,  and  only  if  X  >  0.  Computing  X  from 
Eq.  (66)  yields,  after  simple  manipulations,  the  inequality 


1  T  l-2v 

^3  «(l+v) 


<  0 


(71) 


where  the  upper  or  lower  signs  in  Eqs.  (69)  to  (71)  are  to  be  used  consistently. 
The  limitations  or  >  0  and  V  <  ^  Indicate  that  the  lower  sign  never  leads  to  a 
valid  solution.  Using  the  upper  sign  one  obtains  the  requirement 
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(72) 


a  < 


l~2v 

1+v 


and  the  corresponding  velocity 


(1  H-  SSt/^f 

[1  ^  (i^)l 


(73) 


The  stresses  s^^  and  at  the  front  have  the  ratio 


^i. 

*^1  “  *"2  ”  /3 


(74) 


The  fxinction  f  being  arbitrary,  it  may  be  selected  as  a  step  function, 

Fig.  5b.  The  function 

f(y  -  5t)  «  H(5t  -  y)  (75) 


vanishes  for  positive  values  of  (y  -  ct)  and  is  equal  to  unity  for  negative 
values  of  (y  -  ct).  The  discontinuities  As^  and  in  the  stress  history 


®1  “  ^®1 

(76) 

-  AJ^^  H(ct  -  y) 

(77) 

then  satisfy  Eq.  (74)  provided 

^^1  3er 

.A 

—  »  •/ 

(78) 

The  corresponding  velocity  to  be  obtained  from  Eq.  (63)  is 

V  »  Av  H(ct  -  y) 

(79) 

where 

Av  1  +  2ar  /3 

■  3p3 

(80) 
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Equations  (76)-(80)  give  the  relations  for  an  inelastic  shock  front 
entering  a  stressless  region,  a  case  which  will  be  utilized  in  the  construction 
of  solutions  in  Section  III.  For  completeness  it  is  noted  that  in  general,  the 
region  ahead  of  the  front  need  not  be  stress  free,  but  Eq.  (74)  must  be  satisfied 
ahead  of  the  front. 

The  above  investigation  of  possible  shock  fronts  was  based  on  the  premise 
that  T  *  0.  The  fact  that  no  inelastic  fronts  are  possible  when  t  /  0  is 
demonstrated  in  Appendix  A,  so  that  the  discontinuity  described  in  this  section 
is  the  only  one  which  can  occur  with  inelastic  deformations. 

Summarizing,  it  has  been  demonstrated  that,  for  values  of  a  and  v  which 
satisfy  Eq.  (72),  a  plane  pressure  discontinuity  will  propagate  with  velocity 
c  given  by  Eq.  (73)*  In  order  to  occur  in  the  solution  of  the  steady  state 
problem.  Fig.  4,  the  front  must  be  inclined  at  such  an  angle  9  that  the  horizontal 
component  of  the  velocity  c  equals  the  velocity  V  of  the  load  on  the  surface. 

The  angle  is  obtained  from  Eq.  (73) 
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c.  Regions  emd  shock  fronts  vlthout  inelastic  deformation 


In  a  region  where  no  instantajieous  inelastic  deformation  occurs,  the  strain 
rates  are  defined  by  the  purely  elastic  relations,  Eq.  (11 ),  and  the  stresses  are 
subject  to  the  inequalities 


In  addition,  the  equations  of  motion,  Eqs.  (l8),  (19 )  hold.  To  obtain  the 
differential  equations,  one  could  proceed  in  the  same  manner  as  in  subsection  a. 
However,  it  is  not  necessary  to  do  so,  because  the  resulting  differential  equations 
must  obviously  follow  from  Bqs.  (38)  t)y  making  the  following  two  changes; 

1)  The  last  equation  is  to  be  omitted  because  it  represents  the  yield 
condition  F  «  0,  which  does  not  apply. 

2)  To  account  for  the  change  in  the  stress  strain  law  from  Eq.  (8)  to 
Eq.  (11),  L  0  is  to  be  introduced  into  Eqs.  (38). 
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Tlie  equationsare  linear  and  homogeneous  so  that  the  derivatives  of  the 

t  t  I 

stresses,  Sj^  ,  s^  ,  vanish,  unless  the  determinant  of  Eqs.  (85)  equals  zero. 
In  spite  of  the  fart  that  the  coefficients  in  Eqs,  (85)  contain  y,  the  value 
of  the  determinant  is  independent  of  the  value  of  the  angle  y.  The  determinant 
of  Eqs.  (85)  vanishes  when  X  is  a  root  of 


UX(1-2X)  [1  +  (l-2X)(l-2v)]  -  0 


(86) 


Equation  (86)  has  two  significant  roots. 


and 


1-v 

l-2v 


> 


(87) 

(88) 


and  one,  X  =*  0,  which  may  be  shown  to  be  trivial.  Substituting  the  two  roots 
Xp  and  X„  into  Eq.  (37)  furnishes  two  locations 


where  the  determinant  of  Eqs.  (85)  vanishes,  and  Cp  ,  Cg  are  the  velocities 
of  P-  and  S-waves,  respectively.  In  any  location  cp  /  cpp  or  cpg  the  derivatives 

I  I  I 

Si  ,  Sg  j  vanish,  so  that  the  stresses  must  remain  constant  everywhere, 
except  at  the  locations  (pp  and  jpg  . 

The  angles  cpp  and  ^g  being  the  potential  locations  of  elastic  P  and  S 
shock  fronts,  respectively,  it  is  known  that  discontinuities  in  stresses  and 
velocities  may  occur  at  these  locations  and  may,  therefore,  be  part  of  the 
comple^e  solutions  to  be  obtained  in  Section  III.  The  following  pertinent 
details  for  these  fronts  will  be  required  subsequently. 
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(1)  The  P-front 


Designating  the  changes  in  the  various  quantities  at  the  front  by  the  symbol 
t,  the  discontinuities  in  the  stresses  cr^j  ,  (normal  and  tangential  to  the 

front,  respectively)  and  in  the  component  u^j  of  the  velocity  (normal  to  the  front) 
are  proportional, 


1 


.  V  .  -1 

1-v  •  pCp 


(91) 


No  other  discontinuities  can  occur. 


The  changes  Aa^j  and  Ao^  are  of  course  restricted  by  the  fact  that  the  in¬ 
equalities  (84)  for  the  stresses  must  be  satisfied  on  either  side  of  the  front. 

No  general  study  of  this  restriction  is  required,  but  in  Section  III  it  will  be 
necessary  to  know  if  a  P-front  is  possible  when  the  stresses  and  velocities  ahead 
of  the  front  tpp  vanish.  In  this  case  the  conditions  (84)  are  satisfied  ahead  of 
the  front,  F  »  ■  0.  To  check  behind  the  shock,  it  is  noted  that  the  stresses 

ACj^,  and  Aa^  are  not  only  the  total  stresses,  but  they  are  also  the  principal 
stresses,  Ad^j  »  0^  ,  Ac^  =  dg  .  After  computation  of  and  ,  one  finds  the 
necessary  conditions  for  a  P-front 


d^<0 


(92) 


<x  > 


-  /3 


l-2v 

1+V 


(93) 


A  compressive  shock  front  of  arbitrary  strength  d^^  <  0  in  the  location  cp  »  (pp  is 
therefore  possible  if,  and  only  if,  the  inequality  (93)  on  a  is  satisfied.  The 
angle  v  and  the  quantity  8,  Eq.  (40),  inuaediately  following  the  front  are 


Y  = 


TT 

2 


3  -  3 


(94) 

(95) 


Attention  is  drawn  to  the  fact  that  the  inequalities  (93)  ond  (72)  which 
permit,  respectively,  an  elastic  or  an  inelastic  pressure  discontinuity  to  enter  a 
stress  free  region,  are  mutually  exclusive,  but  complementary.  In  other  words,  for 
any  combination  of  a  and  v, one,  but  only  one,  of  the  two  types  of  fronts  exists. 
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(2)  The  S-front 


At  an  S-front  discontinuities  occur  only  in  the  shear  stress  =«  t 

NT 

and  in  the  tangential  velocity  u^  .  The  changes  are  proportional, 


At  :  Au^ 


(96) 


In  addition,  the  inequalities  (8U)  must  again  be  satisfied  ahead  and  behind 
the  front.  Checking  the  situation  if  the  region  ahead  of  the  front  is  stress 
free,  Eqs.  (84)  are  again  satisfied  ahead  of  the  front.  Behind  the  front  the 
stresses  are 


T  .  AT  ,  .  0  (97) 

The  invariant  Jg  may  be  written 

^2^4*  °N°T 

and,  because  in  this  last  equation  all  stresses,  except  t,  vanish,  the  condition 

P  «  Jg  -  or^J^  <  0  (99) 

is  violated.  Thus  an  S-front  can  not  enter  a  stress -free  region. 

Some  details  about  S-front  locations  where  stresses  ahead  of  the  front  do 
not  veuiish  will  be  required.  Consider  specifically  the  possibility  of  such  a 
front  at  a  point  where  the  equal  sign  in  the  first  condition  (84)  applies  ahead 
of  the  front, 

F  ■  Jg  -  »  0  (100) 

which  indicates  that  the  material  is  at  the  verge  of  inelastic  deformation. 
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Let  the  shear  stresses  just  ahead  of  the  front,  for  be  t,  and 

'  b 

those  behind  the  front,  for  be  T.  The  invariant  and  the  deviators 

with  respect  to  the  N  and  T  directions,  8„  and  s„  ,  respectively,  are  equal  at 

Pi  jL 

9^"^  and  Noting  that  the  equality  (lOO)  is  satisfied  ahead  of  the  front, 

O  b 

it  is  clear  that  the  inequality  Eq.  (99)  requires 

|t|  <  1^1  (101) 

The  largest  possible  change  At  *  T-T  occurs  therefore  when  t  »  -  T,  in  which 
case  Eq.  (100)  is  satisfied  also  for  It  is  useful  to  consider  this  case 

in  terms  of  the  principal  stress  variables  used  in  liubsection  a,  i.e.  using  the 
angle  y  and  the  quantities  ,  s^^  and  |5  as  variables.  Figure  6  shows  the 
direction  of  the  major  deviator  ahead  of  the  front  at  an  angle  y  to  the 

S-front.  The  state  of  stress  ahead  of  the  front,  y>  ^  corresponds  to 

the  values  of  t,  s^j  and  s^  which  apply  in  this  location.  Behind  the  front  the 
state  of  stress  is  defined  by  t  ■  -  t,  Sj^  ,  s^  ,  which  stresses  define  changed 
values  Y;  ^  con5)uting  these  values  by  the  conventional  relations 

it  is  found  that  8^,82,  being  even  functions  of  t,  are  necessarily  equal  to 
81  ,  Sg  ,  respectively,  y,  being  an  odd  function  of  t,  changes 

Y  ■  TT  -  Y  (102) 

Therefore,  a  change  in  shear  from  t  to  t  »  -  T  at  9g  does  not  change  the 
values  of  the  variables  s^  ,  Sg  ,  or  3,  but  only  the  values  of  the  angles 
Y  and  0.  The  latter  becomes 

at  9^"^^  e-9g-Y»9s+Y-fT  (103) 

These  changes  in  y  or  6  occur  if  the  stresses  satisfy  the  equality  (100)  ahead 
of,  and  behind  the  front. 
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Ill  CONSTRUCTION  OF  SOLUTIONS 


In  Section  II  a  nunfiier  of  partial  solutions  were  obtained  from  which  the 
solution  of  the  complete  boundary  value  problem  is  now  to  be  constructed. 

Section  Ila  gives  the  differential  equations  for  the  determination  of  the 
stresses  and  velocities  in  inelastic  regions;  from  Section  lie  it  J  s  known 
that  all  unknowns  in  elastic  regions  are  constants,  except  for  discontinuities 
of  a  prescribed  nature  at  the  locations  fg  and  «pp  ,  In  addition,  there  may  be 
a  shock  front  with  inelastic  deformation  at  a  location 

As  mentioned  in  the  last  two  paragra5>h8  of  the  Introduction,  steady-state 
problems  of  the  type  studied  here  need  not  have  unique  solutions.  However,  it 
may  be  possible  to  eliminate  excess  solutions  by  specifying  that  the  steady-state 
solution  desired  should  be  the  asyn5)totic  solution,  if  any,  of  the  problem  of  an 
expanding  load  (Fig.  7)  applied  on  a  half-space  initially  at  rest.  This  additional 
condition  is  invoked  here  and  furnishes  a  vital  boundary  condition  for  the  solution 
through  the  reasoning  which  follows.  It  is  known  that  the  partial  differential 
equations  of  the  transient  problem,  Fig.  7,  are  hyperbolic  in  elastic  and  in  in¬ 
elastic  regions »  The  characteristic  velocities  U  under  elastic  conditions  are 
U  *  c_  and  U  «  c-  ,  while  those  in  the  inelastic  case  are  functions  of  the  stresses, 
but  subject  to  the  Inequality  U  <  Cp  .  The  hyperbolic  character  of  the  differential 
equations  euid  the  inequality  have  been  demonstrated  by  Mandel  [8J  for  a  general 
class  of  elastic-plastic  materials  governed  by  a  plastic  potential.  This  result 
applies  here.  The  largest  characteristic  velocity  being  Cp  ,  one  can  conclude 
that,  in  the  non-steady -state  stpcrsdsmlc  problem,  V/cp  >  1,  Fig.  7,  all  unknowns 
vanish  ahead  of  a  front  inclined  at  an  angle  <pp  corresponding  to  Cp  ,  such  that 
one  has  a  boundary  condition  for  the  solution  of  the  steady-state  problem 

for  <p<9p:  a^  =  s^  =  u  =  v  =  0  (10^4) 
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Additional  boundary  conditions  apply  at  the  loaded  surface  where  the 
pressure  p^H(Vt  -  x)  is  applied.  At  this  surface  one  of  the  two  principal 
stresses  must  be  vertical  and  equal  to  -  #  so  that  there  are  two  alternative 

boundary  conditions.  Either 


*1 

^1  “  ®1  3  ’^l  "  “  ^o  ^  "  2  (105) 

or 

cfg  =*  Sg  +  ^  Jl  *  -  Po  Y  *  0  or  TT  (106) 

It  can  easily  be  shown  that  the  nature  of  all  equations  in  Section  II  is  such 
that  p^  will  appear  as  an  external  factor  in  the  solution*  for  the  stresseb  and 
velocities^  while  the  nondimensional  quantities  0,  0  and  y  are  independent  of  p^  . 
This  simplification  is  due  to  the  homogeneou*  nature  of  the  plastic  potential; 

Eq.  (4);  and  would  not  apply  if  Eq.  (1),  allcndng  for  cohesion;  is  specified. 
Therefore;  only  the  case  «  1  need  be  considered;  so  that  Eqs.  (105 );  (106) 
become 

<Ti  -  8i  Ji  -  -  1  ,  Y  «  I  (107) 

or 

OgWSg+^Jl*-!  ;■  Y"  (108) 


At  this  point  is  must  be  stressed  that  no  uniqueness  or  existence  theorem 
for  transient  problems  is  available  for  elaatic-plMtic  materials.  Although  the 
boundary  conditi  on  cXl.iuinoLv€  5  certain  excess  solutions  of  the  steady-state 

problem;  Fig.  3,  which  clearly  are  not  eisymptotic  solutions  of  the  transient 
problem;  Fig.  "the  I'emaining  solutions  of  the  steady-state  problem  may  still 
not  be  unique;  because  the  original  transient  problem  may  not  have  a  imique 
solution.  In  constructing  solutions  it  must  be  attempted  to  consider  all  con¬ 
ceivable  possibilities;  but;  in  view  of  the  numerical  procedures  necessary;  it 
is  not  an  absolute  proof  of  uniqueness  if  Just  one  solution  is  actually  found. 
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In  order  to  have  confidence  that  the  solutions  obtained  are  the  physicall;/’ 
meaningful  ones,  even  if  others  should  exist,  the  solutions  will  be  considered 
as  functions  of  the  basic  physical  parameters  v,  a  and  of  v/cp  ,  in  the  expectation 
that  there  should  be  a  continuous  transition  in  character  and  in  the  numerical 
values  at  the  solutions.  The  principle  that  the  character  of  the  solution  should 
change  smoothly  ae  a  function  of  the  parameters  is  also  extremely  helpful  in  the 
formulation  of  the  solutions.  As  a  starting  point  one  can  explore  the  existence 
of  a  range  in  the  above  parameters  where  the  elastic  solution  applies  and  from 
there  continue,  step  by  step,  into  further  ranges. 

Applying  this  gradual  approach,  one  arrives  at  the  conclusion  that  the 
occurrence  of  discontinuous  fronts  at  the  transition  from  the  stressless  to  the 
stressed  state  in  the  elastic  situation  must  also  apply  in  the  more  general  case, 
at  least  for  values  of  the  parameters  close  to  those  where  the  elastic  solutions 
are  valid.  The  first  atten^jt  will  therefore  be  the  construction  of  solutions  with 
an  initial  discontinuity  at  the  arrival  time,  and  the  possibility  of  continuous 
behavior  at  arrival  will  be  considered  subsequently  to  demonstrate  imiqueness, 
and  also  in  situations  where  the  assumption  of  an  initial  discontinuity  does  not 
lead  to  a  solution. 

In  accordance  with  the  above  approach,  one  expects  that  the  properties  of 
the  initial  discontinuity  will  govern  the  character  of  the  solution  as  a  function 
of  the  parameters  v,  or  and  V/c^  .  The  existence  and  the  nature  of  the  initial 
discontinuities  depend  on  v  and  a  only,  such  that  these  two  parameters  play  a  more 
ingjortant  role  than  V/cp  ,  and  a  preliminary  classification  of  the  ranges  can  be 
based  on  v  and  a  only.  In  the  permissible  range  for  these  parameters,  0  <  v  <  ^ 
and  0  <  ®  <  1/  /I?,  there  is  according  to  Section  lie  a  Range  1,  Fig.  9a>  defined 
by 

(109) 


30 


^3 


FIG.9C.  SUBRANGES  FOR  ^Cp  “  2 
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•  numerical  results 

FlG.9d  SUBRANGES  FOR  V/^p  =  1.25 
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where  an  elastic  P-front,  but  no  other  discontinuity  can  enter  a  stress  free 
region,  while  for 


/3  «  < 


l-2v 

1+v 


(no) 


only  a  conqpreBslve  discontinuity  with  Inelastic  deformation  may  enter  a  stress- 
free  region.  The  total  range  where  the  latter  Inequedity  applies  is  subdivided 
in  Ranges  II  and  III,  depending  on  whether  the  velocity  c  of  the  inelastic  front 
is  A.arger  or  smaller  ichan  the  velocity  Cg  of  elastic  shear  waves,  respectively. 
Range  III,  where  c  <  Cg  ,  applies  if 


/3  or  <  -  2  +  -pr:?-==: 

■/T(l+^ 


(111) 


while  Range  II,  c  >  Cg  ,  applies  if  Eq.  (ni)  is  violated.  The  reason  for  this 
division  will  be  seen  later. 


a*  Range  la 

As  indicated  above,  the  first  step  in  the  construction  of  solution  is  the 
determination  of  the  range,  designated  Range  la,  in  which  entirely  elastic 
solutions  exist.  In  such  a  solution  a  P-wave  enters  a  stressless  region  and  an 
entirely  elastic  solution  is  possible,  if  at  all,  only  in  a  range  entirely  within 
Range  I,  Fig.  9^.  The  stresses  in  an  eleistic  half -space  due  to  a  superseismically 
traveling  uniform  surface  pressure  are  given  in  Appendix  B.  There  is  a  uniform 
state  of  stress  between  the  P-front  and  S-front,  and  again  a  uniform,  but  different 
state  of  stress  between  the  S-front  and  the  surface.  Fig.  10.  The  two  uniform 
states  of  stress  must  satisfy  the  inequalities  (84) 

^  <  0  *(112) 
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DIRECTION  OF  0? 


S-FRONT 


F  I  G  .  10  CONFIGURATION  FOR  RANGE  IQ 


FIG.  I  I  CONFIGURATION  FOR  RANGE  Ib 
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There  is  no  need  to  check  the  validity  of  these  inequalities  in  the  region 
between  the  P  and  S  fronts  because  this  has  already  been  done  in  Section  Tic 
where  the  existence  for  the  P-front  in  Range  I  was  proved  However,  Eqs.  (112) 
must  be  considered  for  rf  >  tp  >  9  .  At  the  S-front  a  state  of  pure  shear  is 
added  to  the  state  of  stress  for  9  <  cpg  .  This  can  not  change  the  first  in¬ 
variant  ,  so  that  only  the  condition  F  <  0  requires  checking. 


Substitution  of  Eqs.  (179 (I80)  from  Appendix  B  into  this  condition 


results  in  the  inequality 


,2  .  1  A-2vJ  .  l' 


a  > 


-If  VrnJ 


cos  29g 


(113) 


The  inequality  (113)  defines  Range  la  where  the  response  is  entirely  elastic. 

The  range  is  a  function  of  Poisson's  ratio  and  of  V/cp  >  1,  and  its  boundary 
can  be  found  by  using  the  equal  sign  in  Eq.  (113)*  Figures  9t>“d  show  that  these 
boundaries  end  at  ..ne  one  between  the  principal  Regions  I  and  II,  the  endpoint 
being  defined  by  the  relation 


Figures  Jb-d  show  Range  la  covering  ne^ly  all  of  Range  I,  while  in  Fig.  ye 

Range  la  actually  covers  all  of  Range  I.  Using  .ilqs.  (113 )»  (H^)  the 

limiting  value  a  =  one  finds  the  critical  value  V/Cp  =  I.06I  below  which 

V12  ^ 

Range  la  covers  all  of  Range  I.  The  stresses  in  Range  la  are  entirely  elastic 
and  are  given  by  the  simple  relations  listed  in  Appendix  B. 


"b.  Range  Ib 


It  was  found  above  that  entirely  elastic  solutions  exist  only  when  tlie 
inequality  (113)  is  satisfied.  The  remainder  of  Range  I,  i.e.  the  range 


®  2{?3  -  <(>p) 

2 

cos  acpg 


+ 


1 

3 


(117) 


will  be  designated  as  Range  Ib.  In  this  range  the  solution  can  no  longer  be 
entirely  elastic  and  must  therefore  contain  at  least  one  location  with  inelastic 
deformation. 


Using  the  expected  continuity  of  the  character  of  the  solutions  as  a  guide, 
the  solution  in  this  range  tc  start  agaj.n  with  a  discontinuity  which, 

according  to  Section  lie,  can  only  be  an  elastic  P-front  located  at  fp  .  Using 

/  x  (+) 

Eq.  (91)  for  the  stress  changes  at  the  front,  one  finds  that  for  '  9  “  9p  the 

l*2v 

inequality  F  <  0  is  satisfied,  provided  the  special  case  /3  or  ■  is  excluded 
for  separate  consideration.  Having  recognized  that  the  solution  must  contain  an 
inelastic  region,  where  F  ■  0,  a  further  elastic  stress  change  must  occur,  which 
is  possible  only  at  the  S-front.  The  appropriate  change  in  the  state  of  stress 
at  the  S-front  has  been  obtained  in  Appendix  C,  in  terms  of  the  as  yet  unknown 
stress  discontinuity  Ao  at  »p  . 

For  <  9  <  9s” ^  ' 


«  A<J  ,  0  «  3  ^  Y  "  I  ^  ®  "  i 


while  for  9 


1+v 

1-v 


ha 


(118) 


(119) 

(120) 


( *4" } 

The  symbol  (  +  )  in  <pp  '  indicates  a  value  infinitesimally  larger  than  cpj, 
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(121) 


8 


1 


11-2V)(94-.1) 

5(l-v) 


ha 


Y  =  I  +  |a| 

®  “  i  ±  Ul 


The  quantity  fi  is  obtained  from 

cos  2fi  -  ^  cos  2(«pp  -  qjg) 
and  is  subject  to  the  inequality 


(122) 

(123) 


(124) 


9s  >  > 'Ps  “  (^25) 

r  1-2V 

The  special  case  / 3  »  »  1+v”  discussed.  In  this  case  F  =  0 

(+) 

is  satisfied  already  for  cp  «  «^p  so  that  the  possibility  of  an  inelastic 

region  no  longer  requires  a  shear  front  at  cpg  .  However,  a  change  in  shear 

leadi  g  again  to  a  state  with  F  =  0  is  still  possible.  Both  possibilities  are 

actually  included  in  Eqs.  (119)“(124).  The  special  case  simply  means  that  one 

of  the  two  values  ®(9s'*’^)  is  equal  to  9(cpp^^)  given  by  Eq.  (Il8).  The  fact 

that  in  the  special  case  an  inelastic  region  may,  in  principle,  occur  in  the 

(  + ) 

range  cpp  '  <  cp  <  cpg  is  of  no  consequence,  because  such  a  solution  would  violate 
the  necessary  continuity  of  the  configurations. 


The  resxilts  obtained  so  far,  and  further  steps  required,  are  best  discussed 

in  terms  of  the  angle  9  in  various  locations,  illustrated  in  Fig.  11.  The 

direction  of  the  principal  stress  between  the  P  and  S-fronts  according  to  Eq.  (Il8) 

{  +  ) 

is  normal  to  the  P-front,  while  for  cp  >  9g  Eq.  (123)  defines  9.  Because  of 
the  inequality  (125 )^  ®(9g^^)  is  less  than  ^  but  more  than  0,  regardless  of  the 
sign  of  6.  According  to  Section  II,  there  is  no  further  possibility  for  a  change 
in  9  as  required  to  arrive  at  the  surface  value  0(n)  =  ^  ;  (0  or  n),  except  one 
or  more  inelastic  regions  for  9  >  cp„  . 
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Using  the  values  of  0,  ,  s^  ,  Y  6  defined  by  Eqs.  (119)-(123}  the 

results  of  Section  Ila  are  now  to  be  used  to  find  and  determine  the  history  of 
the  stresses  and  particularly  of  the  angle  0.  If  a  region  can  be  found  during 
the  forward  integration  where  either  of  the  values  B  =  ^  (or  0,  or  tt)  is  obtained, 
the  integration  is  terminated.  From  the  point  of  termination  to  the  surface  an 
elastic  region  of  no  change  is  selected  such  that  the  surface  condition  for  0  is 
then  satisfied.  During  this  integration  the  unknown  value  Aa  in  Eqs.  (120),  (121) 
is  a  common  factor  in  all  stresses,  so  that  the  integration  will  give  a  principal 
stress  at  the  surface,  which  contains  this  factor,  which  is  selected  to  satisfy 
the  boundary  condition,  Eq.  (107  or  108),  0^  ^  -  1. 

The  use  of  the  solutions  derived  in  Section  Ila  for  inelastic  regions  is 

( + ) 

quite  straightforward.  From  the  values  of  0,  y  tpg  f  potentieil  starting 
points  <p^  of  inelastic  regions  are  located  as  roots  of  the  determinantal  equation 
(40).  Next  it  must  be  verified  that  GL/V,  Eq.  (53);  is  negative,  If  this  is  so, 
Eqs.  (44)-(47)  are  used  to  determine  the  solution  by  forward  integration,  con¬ 
tinuously  checking  the  sign  of  GL/v.  The  integration  can  be  continued  until  GL/V 
changes  sign,  but  may  be  stopped  at  any  desired  location  epg  .  When  an  angle 
0  =  ^  (or  0,  or  tt)  for  the  direction  of  the  principal  stress  is  obtained,  a 
solution  to  the  problem  has  been  found. 

The  configuration  considered  was  successful  and  led  to  just  one  solution  of 
the  problem.  The  upper  sign  in  Eq.  (123)  and  the  case  0  =  ^  furnished  the  solution, 
but  it  is  suspected  that  this  may  not  be  so  when  V/Cp  >  1  is  close  to  unity.  The 
matter  of  possible  alternative  configurations  which  might  lead  to  solutions  is 
discussed  later  in  this  section. 

It  is  noted  that  Range  Ib,  which  does  not  occur  at  all  if  v/cp  <  I.061, 
applies  even  for  other  values  of  V/cp  only  in  a  minute  portion  of  the  overall 
range  of  v  and  a,  as  can  be  seen  from  Figs.  (9b-d). 
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Accordi.ng  to  the  definition  of  ranges  at  the  beginning  of  this  section, 
an  initial  inelastic  discontinuity,  but  no  other,  is  possible  in  Range  II. 
J'\irther,  in  this  range,  the  location  9  of  this  discontinuity,  defined  by 
Eq.  (81),  is  such  that 

9  <  9g  (126) 

Range  II,  vhich  is  the  one  of  major  interest,  is  defined  by  the  inequalities 

(110,  111), 

— L=r-2</3a<^  (127) 

/2(l+vy 

In  that  portion  of  Range  II  which  adjoins  Reinge  la.  Figs.  9^-e,  one  expects 
that  the  solutions  after  starting  with  an  inelastic  front  of  discontinuity  at  9 
will  remain  entirely  elastic.  The  range  in  which  such  solutions  apply  and  the 
values  of  the  stresses  are  obtained  in  Appendix  D.  This  range  is  designated 
Range  Ila,  and  the  stresses  are  found  in  closed  form,  the  configuration  being 
shown  in  Fig.  12. 

The  discontinuity  to  in  the  normal  stress  at  the  front  is 

-  cos  29„ 

ta  -  - : - - 5 -  (128) 

(l-R)  cos  (9  -  9g)  +  (1+R)  cos  «Pg  -  1 

where 

R  -  ~  (129) 

1  +  2ar  /3 


U6 


FIG.  12  CONFI  GURATION  FOR  RANGE  ILQ 


FIG. 13  CONFIGURATION  FOR  RANGE  H b 


is  satisfied.  The  boundary  separating  Region  Ila  from  the  remainder  of  Region  II, 
designated  Region  Ilb,  is  found  by  using  the  equal  sign  in  the  above  relation. 
Figux'es  9^1  to  e  show  typical  curves  for  some  values  of  V/cp  .  These  figures  in¬ 
dicate  that  Range  Ila  covers  only  a  quite  small  portion  of  Range  TI,  except  in 
the  atypical  case  when  V/cp  is  only  slightly  larger  than  unity,  Fig.  9e. 
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d.  Range  lib 


In  Range  II ^  but  outside  Range  Ila,  the  solutions  are  expected  to  stai't 

with  an  inelastic  pressure  front  at  but  additional  inelastic  regions  must 

now  occur.  In  the  vicinity  of  the  boundary  towards  Region  Ib,  continuity 

requires  sirailai’  configurations,  as  shown  in  Fig.  13*  Behind  the  inelastic 

front  the  stresses  will  be  uniform  with  a  shear  front  at  (pg  ,  and  a  region  of 

inelastic  deformation  in  a  location  «pg  <  ^  <  n.  The  discontinuity  in  shear  At 

(  +  ) 

must  be  such  that  the  yield  condition  F  »  0  is  satisfied  for  9  >  «Pg  • 

However,  at  points  remote  from  the  boundary  between  Regions  Ib  and  Ilb 
alternative  configuratic is  might  occur  and  must  be  considered  as  possibilities 
in  the  numerical  analysis.  In  the  configuration  shown  in  Fig.  13  the  possibility 
At  »  0  could  furnish  a  solution,  or  inelastic  regions  may  exist  in  locations 
9  <  «p  <  as  shown  in  the  alternative  Figs,  ik  and  15>  where  shear  discon¬ 

tinuities  At  ^  0,  may  occur,  or  not,  At  ■  0.  Further,  the  configuration  shown 
in  Fig.  15  may  have  a  subrange  where  the  discontinuity  At  is  such  that  elastic 
conditions  F  <  0  are  created  and,  therefore,  constant  stresses  occur  for  tp  >  • 

Disregarding,  for  later  discussion,  solutions  without  initial  discontinuity, 
but  allowing  inelastic  regions  to  split,  this  exhausts  all  possibilities  to  be 
studied.  The  numerical  analysis  by  con^juter  furnished  only  solutions  having 
the  configuration  of  Fig.  I3.  The  search  for  roots  of  the  determlnantal  equation, 
giving  starting  points  of  inelastic  regions  never  furnished  a  root  for  «p  <  <pg  . 

The  following  statement  summarizes  the  situation.  The  initial  change  from 
vanishing  to  nonvanishing  stresses  occurs  at  an  inelastic  front  with  an  as  yet 
undetermined  compressive  discontinuity  A<y  in  the  principal  stress  <?^  .  This  front 
is  followed  by  a  region  of  constant  stress. 
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FIG.  I  5  alternative  configuration 


<72^  ■  Lc 
P  -  3 


for  ^ 
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For  locations  9  >  ?pg  there  ore  two  alternatives.  If  no  discontinuity  in 
shear  occurs,  At  «  0,  Eqs.  (13^)  apply  also  for  <p  «  while  the  angle  y  is 

V(q?^‘^b  =•  <Pg  -  9  +  I  (135) 

However,  if  a  shear  discontinuity.  At  /  0,  occui's,  Eqs.  (102),  (I03)  give,  using 
Y  «  cpg  -  9  +  I  , 
for  : 

■  Ao 

0*3 

tr  («6) 

Y  -g 

e  .  2,g  -  ^  -  2 

Equations  (I36)  and  the  alternative  values  for  At  ■  0  are  the  starting  points 
for  numerical  integrations  which  are  tc  be  carried  out  in  the  manner  described 
for  Range  lb. 


e.  Search  for  inelastic  solutions  without  initial  discontinuity 

In  Ranges  I  and  II  solutions  were  constructed  where  the  initial  change, 
from  vanishing  to  nonvanishing  stresses,  occured  as  a  shock,  either  elastic  at 
cpp  ,  or  inelastic  at  9.  While  the  principle  of  the  continuity  of  solutions 
makes  the  solutions  obtained  plausible,  it  is  desirable  to  investigate  if 
solutions  which  start  smoothly  exist. 
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V 


The  differential  equations  in  elastic  regions  permit  definitely  no  smooth 
change  in  stress  fo’^  a:?3erseisnic  velocities  V,  such  that  only  the  inelastic  case 
is  considered. 


If  a  smooth  inelastic  solution  starting  from  vanishing  stresses  in  a  location 

cp^  exists,  an  asymptotic  study  of  the  appropriate  differential  equations  for 

a ,  =  8.  =•  J,  -+  0  in  the  vicinity  of  cp  must  describe  this  solution.  In  order  to 
J  J  1  o 

be  physically  sensible,  the  angle  y  in  the  vicinity  of  cp^  must  be  well  behaved 
and  may  be  considered  a  constant  in  the  range  <  9  <  <Pq  +  e  where  e  is  small. 

The  quantity  jL/v  <  0  must  not  vanish,  otherwise  the  region  is  not  inelastic  as 
postulated.  There  are,  however,  two  possibilities  for  the  behavior  of  GL/v. 

In  the  limit  cp  -♦  cp^  ,  the  function  GL/v  may  be  finite  and  well  behaved,  in  which 
case  it  may  be  considered  a  constant  near  <p^  j  alternatively,  GL/v  may,  in  the 
limit,  be  infinite. 


The  first  possibility,  where  GL/v  in  the  limit  may  be  replaced  by  a  constant 

is  easily  proved  to  be  impossible.  Following  the  previous  reasoning  in  Section  II, 

solutions  in  an  inelastic  region  exist  only  if  the  determinant  of  Eqs.  (38) 

vanishes,  in  which  case  four  of  the  five  unknowns  will  depend  on  the  fifth.  In 

the  limit  -*  0  the  last  Eq.  (38)  becomes  trivial,  0  ■  0.  The  last  terms 

of  the  other  equations  vanish,  because  GL/V  is  finite  and  products  of  GL/V  and 

Sj  or  in  the  limit  Eire  therefore  zero.  The  remaining  four  equations  are  then 

identical  with  the  four  Eqs.  (85)  in  the  elastic  case.  They  have  nonvanishing 

solutions  only  when  <p^  «  cpg  or  cp^  -  cpp  .  However,  the  yield  condition,  F  =  0, 

represented  by  the  last  Eq.  (38),  which  became  trivial,  may  now  not  be  satisfied 

and  must  be  checked.  In  the  vicinity  of  s^  ,  -*  0,  the  ratio  of  these  stresses 

must  obviously  be  the  same  as  at  the  P  or  S  front,  Eqs.  (9I)  and  (98),  respectively, 

obtained  from  the  same  equations.  Based  on  the  discussion  of  the  P  and  S  fronts, 

one  finds  easily  that  the  requirement  F  =*  0  is  not  satisfied,  except  in  the  special 

case  when  the  values  O'  and  v  are  exactly  on  the  boundary  between  regions  I  and  II, 
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where  9  a  qjp  .  However,  in  this  case  one  finds  GL/V  «  0,  and  no  inelastic 
solutions  whatsoever  are  therefore  possible  when  GL/V  at  9^  is  finite. 

The  case  where  |gL/v(  •  as  <p  -♦  9^  remains  to  be  discussed.  The  first 
question  concerns  the  possibility  of  |GL/vj  -*•  •  and  conditions  for  the  occurrence 
of  such  a  singularity.  If  such  a  point  exists  for  some  values  of  y  and  of  the 
ratios  of  and  ,  when  the  latter  are  small,  0,  then  |GL/vj  -♦  •  would  also 
occur  for  the  same  ratios  if  s^  and  are  finite.  Equations  (44),  (45)  and  (47) 
which  applj’’,  would  then  give  infinite  values  for  one  or  more  of  the  derivatives 

f  f 

Sj  ,  .  The  possibility  Igl/v|  -♦  •  exists  therefore  only  in  locations  where  an 

inelastic  front  of  discontinvilty  may  occur  and  the  conditions  required  are  those 
for  such  a  front.  Using  the  results  of  Section  Ilb  for  Inelastic  discontinuities, 
no  smooth  solution  can  exist  in  Region  I,  there  being  no  inelastic  shock  front. 

In  Regions  II  and  III  where  such  a  front  is  possible  at  9,  a  solution  of  the  type 
sought  may  exist, starting  at  9^  »  9  ;  the  necessary  initial  values  of  the  state 
of  stress  being  again  defined  by 


P{«Pq)  *  3 

■  I 


(137) 


If  an  inelastic  solution  in  the  region  9  >  9  actually  exists,  the  determi- 
naintal  equation  (40)  must  be  satisfied  for  9  >  9.  This  is  necessary  because  the 
previous  reasoning  only  implies  that  this  equation  is  satisfied  for  9  »  9.  To 
explore  this  point  an  asymptotic  expression  for  Eq.  (4o)  is  obtained  by  substi¬ 
tuting 


9  > 

e  - 

Y  = 


9  +  « 

3+4 
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(138) 


where  A  and  T]  ore  sraall  quantities.  Retaining  the  lowest  order  terras  in 
the  new  variables  one  obtains  the  expression 

bgtl®  -  4^  .  bj.  (139) 

where  the  quantities  bg  and  b^  are  functions  of  v,  ot  and  V/Cp  ,  given  in 
Appendix  E.  The  quantity  b^  is  always  positive,  while  bg  may  be  positive  or 
negative,  changing  the  character  of  the  equation  radicedly. 

In  Range  II,  i.e.  when  the  inequality  (127)  applies,  bg  is  negative  so 
that  the  equation  has  real  roots  only  for  negative  t.  While  an  inelastic  region 
can  exist  for  ^  <  9  ending  at  ^  with  vanishing  stresses,  no  such  regions  can  exist 
for  9  >  9,  i.e.  in  the  location  of  interest  here.  The  solutions  in  Ranges  I  and  II 
with  an  initial  discontinuity  previously  obtained  are  therefore  unique. 

In  Range  III,  where  the  inequality  (111)  applies,  one  finds  bg  >  0,  so  that 
the  deterrainantal  equation  has  real  roots  for  i  >  0  as  necessary  for  solutions 
without  initial  discontinuity  in  stress.  The  fined,  condition,  GL/V  <  0,  is  also 
satisfied,  because  the  stress  ratios  in  this  region  are  initially  equal  to  those 
for  the  inelastic  front,  where  GL/v  <  0.  All  requirements  are  therefore  satisfied 
and  it  is  concluded  that  in  Range  III,  and  only  in  this  range,  an  inelastic 
solution  without  stress  discontlnviity  exists.  The  details  of  its  determination 
are  given  in  Appendix  E. 

Range  III 

According  to  the  definition  of  ranges,  an  inelastic  shock  front  in  the 
location  9  is  possible,  and  one  can  atten5)t  to  construct  a  solution  starting 
with  this  discontinuity  in  analogy  to  Range  II.  However,  the  computational  search 
for  inelastic  regions,  for  (p  >  9,  was  unsuccessful,  and  the  boundary  conditions  on 
the  surface  can  not  be  satisfied  in  this  fashion.  The  deterrainantal  equation  (1^0) 


is  nonlinear  and  too  complex  to  prove  the  nonexistence  of  roots  in  general. 


However,  the  approximate  Eq.  (139)  furnishes  a  partial  proof,  as  there  is 
obviously  no  root  ii  »  T]  «  0  for  e  >  0. 

The  impossibility  of  finding  a  solution  with  an  initial  discontinuity  is, 
however,  very  satisfactory  because  the  previous  subsection  and  Appendix  E 
indicate  that  in  this  range  a  solution  exists  which  starts  at  (p  ■  without 
discontinuity.  Because  of  their  singular  character  the  differential  equations 
at  and  near  the  stai'ting  point  can  not  be  solved  by  the  numerical  procedure  used 
in  the  other  ranges.  Therefore,  the  asymptotic  solution  obtained  in  Appendix  E 
must  be  applied  for  a  small  range  until  the  solutions  are  sufficiently 

well  behaved  to  return  to  the  numerical  integration  of  the  differential  equations 
in  Section  Ila. 

It  may  be  surprising,  but  beyond  the  fact  that  a  smooth  solution  may  start 
at  (p  «  only  a  qualitative  statement  on  the  asymptotic  solution  is  actually 
needed,  the  numerical  coefficients  derived  in  Appendix  E  need  not  be  used.  It 
is  sufficient  to  know,  from  Eq.  (2l8),  that  in  the  proximity  of  9  «  9  the  value 
of  0  becomes  approximately 


P  -  3  +  A  ^.3  {1^0) 

while 

Y  -  I  4-  n  ,  (I4l) 

where  the  small  quantity  |  Tl|  is  inherently  larger  than  the  neglected  value  |a|  . 
Equation  (219),  and  a  similar  expressivu  for  the  principal  stress  ,  contains 
an  arbitrary  constant  C  •  Thus,  at  a  point,  ®  =*  9  +  e  ,  still  to  be  selected 

O  ’O'O 

as  end  of  the  asyii5)totic  region,  the  veLLue  of  the  principal  stress  ^^^(9^) 
used  as  the  arbitrary  constant  instead  of  .  Choosing  a  value  Tl  =*  Hq  ^  small, 
yet  large  enough  for  the  numerical  integrations  to  work  thereafter,  one  searches 
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for  the  corresponding  value  where  the  detemiinantal  equation  is  satisfied 
by  the  combination  of  3  «  3^  Vq  \  ^  •  The  principal  stress 

at  this  point  can  be  made  equal  to  vmity.  From  this  point  on  integration  p.r-.v.eedB 
exactly  as  in  the  other  ranges.  Due  to  the  fact  that  Eq.  (217’)  defining  1)  has  a 
+  sign,  it  is  necessary  to  include  the  two  possibilities  +  Tl^  . 

The  procedure  or  Mined  weis  found  to  be  successful,  one  and  only  one,  of  the 
integrations  for  +  furnishing  a  solution.  The  stresses  in  the  interval  ^  to 
cp^  Increase  as  (cp  -  To  obtain  their  distribution  the  exponent  n  can  be 

obtained  from  Eq.  (220).  It  is  a  very  smaLl  positive  number,  of  the  order  of 
l/lOO.  The  configuration  of  solutions  in  Range  III  is  shown  in  Fig.  l6. 

The  occurrence  of  solutions  with  and  without  initial  discontinuity  in  stress, 
does  not  break  the  continuity  in  the  character  of  the  solutions.  Even  for  the 
continuous  solutions  the  derivative  of  the  stresses  at  9  is  infinite,  as  at  a 
discontinuous  front,  and  the  numerical  results  indicate  that  the  change  in  stress 
in  the  a8yn5)totic  region  due  to  the  small  exponent  n  is  so  rapid,  that  this 
region  is  practically  indistinguishable  from  a  discontinuity,  see  Fig.  17* 

t 

g.  Simplified  determination  of  velocities  and-  accelerations 

The  basic  relations  in  Section  II  permit  the  numerical  determination  of 
stresses  and  velocities  or  accelerations.  The  integration  for  the  stresses  must 
be  actually  carried  out  to  obtain  the  appropriate  open  consteint  from  the  boundary 
condition  on  the  surface.  The  parallel  integrations  to  find  velocities  and 
accelerations  may  be  avoided,  by  using  the  following  relations,  some  of  which  are 
exact,  while  others  are  only  good  approximations. 

At  all  fronts  of  discontinuity  the  accelerations  are  of  course  infinite, 
but  the  changes  in  velocity  are  given  -  exactly  -  in  terms  of  the  respective 
stress  dlncontinuities,  if  any. 
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FIG.  17  TYPICAL  RESULT  IN  RANGE  g 
( V  =  0  ,  a=0.05  ,  V/Cp=«  2) 


at  9  -  jpp 

•*’^1  - 

at  »  «pg 

at  cp  >  9 

(142) 

(143) 

(144) 


where  the  subscripts  Indicate  normal  and  tangential  directions,  respectively, 
and  the  velocities,  Cj  ,  are  given  by 

c  -Vsin^  (145) 


The  values  of  at  the  pressure  fronts  can  be  taken  directly  from  the 
numerical  computations  for  the  stresses.  The  value  AT  can  easily  be  computed 
from  the  values  and  directions  of  the  principal  stresses  ,  for  and 

Of  major  interest  is  Range  Ilb  where 

|At|  -  sin  2(fg  -  ^)  <yi(q»g)  (146) 


In  continuous  elastic  regions  velocities  do  not  change,  while  accelerations 
vanish.  Inelastic  regions  being  very  narrow,  one  may  disregard  tangential 
accelerations  and  changes  in  velocity,  while  the  normal  acceleration  may  be 
assumed  to  be  constant  in  the  region,  giving  a  linear  change  in  velocity.  The 


total  change  in  velocity  AUj|  in  an  inelastic  region  of  extent  Aflp,  may  be  found 
from  the  change,  ,  in  the  principal  stress  at  both  ends  of  the  region 


while  the  acceleration  is 


(146) 
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The  changes  in  velocities  at  the  front,  given  exactly  by  Eqs.  (l42)-(lU4), 
ore  always  much  larger  than  those  in  regions  given  by  Eq.  (14? ),  such  that  the 
simplification  of  using  Eqs.  (147),  (l48)  is  quite  satisfactory  when  determining 
shock  factors. 

h.  Nuroarical  analysis 

In  Ranges  Ib,  Ilb  and  III  a  numerical  search  for  inelastic  regions,  and 
subsequent  numerical  quadratures  are  required.  In  Section  II  the  basic  equations 
have  been  written  in  a  very  abbreviated  form,  somewhat  concealing  the  con5)lexity 
of  these  relations.  The  solution  of  these  equations  by  hand  computation  would 
be  nearly  impractical,  and  the  computations  were  made  on  an  IBM  7090.  A  common 
program  was  devised,  allowing  for  the  different  configurations  which  may  occur. 

The  inelastic  regions  are  always  quite  narrow  as  functions  of  only  a  few 
degrees,  and  become  even  narrower  as  v/cp  becomes  large.  It  was  therefore 
necessary  to  vary  the  intervals  of  (p  in  the  search  and  in  the  quadratures.  For 
V/Cp  <  2  intervals  of  I/5OO  rad.  were  used,  while  for  V/Cp  ■  5  intervals  of 
1/10000  ):ad.  were  selected. 

The  results  obtained  are  discussed  in  Section  IV. 
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IV  NUMERICAL  RESULTS  AND  CONCLUSION 


a.  Results 

The  effects  of  a  superseismlcally  progressing  step  pressure  on  the  surface 
of  a  half -space  have  been  obtained  for  an  elastic -plastic  medium  subject  to  the 
yield  condition  (4),  representing  an  inelastic  material  governed  by  internal 
Coulomb  friction.  The  solutions  depend  on  the  elastic  material  parameters  E 
and  and  on  the  additional  parameter  a  in  Eq.  (4).  ff  is  related  to 

the  ajigle  #  of  internal  friction,  using  Eq.  (lO)  of  Ref.  [6], 

sin  f  «  - - ^ - ;  (149) 

In  spite  of  the  lack  of  a  general  uniqueness  and  existence  theorem,  a  unique 
solution  was  obtained  for  each  conibination  of  material  parameters,  sirrface  load 
p,  and  velocity  V/cp  >  1.  There  are,  however,  radically  different  configurations 
depending  on  the  values  of  the  nondimensional  parameters  v,  er  and  V/Cp  .  The 
ranges  in  which  the  various  configurations  apply  haire  been  designated  by  I,  II 
and  III,  where  Ranges  I  and  II  have  been  subdivided  into  Subranges  a  and  b.  The 
values  of  the  parameters  v  and  cr  alone  determine  which  of  the  Ranges  I,  II  or  III 
applies  in  a  particular  case,  as  shown  in  Fig.  9a,  while  the  subdivision  into 
a  or  b  depends  on  the  value  of  V/cp  ,  typical  cases  being  shown  in  Figs.  9b-e. 
These  figures  show  that  Ranges  la  and  lib  cover  most  of  the  total  range  in  v 
and  or,  the  other  ranges  being  of  very  limited  applicability.  Range  la  gives 
entirely  elastic  solutions,  known  from  Ref.  [1],  and  is  not  further  considered. 
The  solutions  found  for  Range  Ilb  are  those  of  prime  interest,  the  other  ranges 
are  somewhat  academic. 

The  parameter  ot  is  inherently  restricted,  or  <  — ■  ,  but  sensible  values 

y  12 

for  the  angle  i  of  internal  friction,  Eq.  (149),  permit  a  further  limitation 
to  the  range  0.10  <  a  <  0.20.  Numerical  resxilts  were  therefore  obtained,  as 
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indicated  in  Figs.  9b-d,  for  combinations  of  V/cp  »  1.25,  ^  ^  ^  J  ^ 

and  a  »  0.10,  0.15,  0.20.  The  values  a  selected  cover  the  range  sin  I  »  0.3  to 
0.7*  Except  for  two  points,  which  feill  in  Range  Ila,  all  these  combinations  are 
in  Range  Ilb.  For  con5>leteneBS  the  result  for  one  case  in  Range  III,  v  ■  0, 
a  m  0.05,  V/Cp  *  2  was  also  obtained. 

Figure  l8  shows  a  typical  variation  of  the  principal  stress  o,  and  of  the 
angle  9  in  the  major  Range  lib.  There  is  a  discontinuous  rise  in  the  principeQ. 
stress  at  the  inelastic  front,  followed  by  a  discontinuity  in  direction  6, 
but  not  in  magnitude  of  ,  at  the  S-front.  There  is  further  a  minor  increase 

in  in  the  inelastic  region  combined  with  a  change  in  direction,  6.  For  unit 

\ 

step  pressures,  p^  *  1,  Table  I  gives  the  values  of  the  principal  stresses  , 

,  and  of  the  angle  9,  and  the  locations  of  the  fronts  for  all  cases 
considered,  which  fall  into  Range  lib. 

Figures  19a,  b  show  and  0  for  the  two  cases,  v  »  l/3,  or  ■  0.10,  V/cp  =  1.25 
and  2.0,  which  fall  into  Range  Ila.  In  these  cases  the  initial  stress  rise  is 
again  at  the  inelastic  front,  9  ■  ^.  There  is  a  change  in  and  0  at  the 

S-front,  but  there  are  no  further  inelastic  regions,  and  no  further  changes  in 
or  0.  The  solution  in  this  range  does  not  require  numerical  integrations,  but 
is  entirely  in  closed  form. 

Range  III  is  only  of  theoretical  Interest,  because  it  applies  only  for 
sin  #  <  0.21,  but  a  typical  case  is  shown  In  Fig.  I7.  There  is  no  discontinuous 
front,  the  solution  starts  smoothly  at  f  but  the  principetl  stress  has  a 
vertical  tangent  and  rises  extremely  rapidly,  nearly  like  a  discontinuous  front. 

b.  Conclusion 

From  a  practical  point  of  view,  tr.e  most  important  conclusion  obtained  by  the 
present  analysis  concerns  the  general  character  of  the  solutions.  The  ntimerical 
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INELASTIC  PRESSURE  FRONT 


FIG. 18  TYPICAL  RESULT  IN  RANGE  Eb 
(  V  =  0,  a  =  O.I5  ,  V/C^a  2) 


67 


TABLE  I.  RESULTS  IN  RANGE  Ilb 


V 

* 

'Ps 

0 

145.55 

1/8 

148.4? 

1/4 

152.49 

0 

159.30 

1/8 

160.89 

1/4 

163.22 

0 

171.87 

1/8 

172.48 

1/4 

173.37 

1/3 

174.26 

Location 


.10  131.53 

.15  128.42 
.20  127.06 


.10  157.82 
.15  155.66 
.20  153.97 


.10  152.10 
.15  150.68 
.20  150.08 


0.10 

0.15 

0.20 


.10 

.15 

.20 


0.10 

0.15 

0.20 


0.10  168.61 


57-7 

56.5 

55.5 


58.26 

57.29 


64.01 

63.17 

62.53 


64 

63 

62 


.97 
.62 
164.38 


2.5 

2.3 


172. 7O 
173.69 
172.64 


173.50 

173.47 

173.45 


174.33 


0.5804 

0.4424 

0.3346 


0.5973 

0.4733 

0.3753 


0.5960 

0.4615 

0.3547 


0.5979 

0.466? 

0.3622 


0.6053 

0.4776 

0.3760 


0.6105 

0.4823 

0.3805 


0.6342 

0.5155 

0.4206 


0.6319 

0.5098 

0.4118 


0.6233 

0.4969 

0.3965 


0.6177 

0.4920 

0.3918 


.6171 

.4908 

.3903 


B  and  C 


O.Oill 

0.7950 

0.7881 


69.32  0.8987  0.5519  47.51 
73.37  0.8791  0.4282  43.47 

76.32  0.8794  C.3395  40.52 


83.46  0.9641  0.5920  41.53 
86.57  0.9743  0.4746  38.42 
87.92  0.9813  0.3789  37.06 


70.77  0.8954  0.5498  67.82 
72.91  0.8943  0.4356  65.68 
74.62  0.8955  0.3457  63.97 


4l  0.5797  65.03 
10  0.4584  63.03 
09  0.3633  61.63 


84.34  0.9791 
85.76  0.9799 
86.36  0.9801 


82.43  0.9703 
O3.22  0.9761 
83.85  0.9789 


.9887 

.9889 

.9893 


87.53  0.9 
88.03  0.9 
88.25 
8 '.91 1 0.9999 


.5797 

.4584 

.3633 


.6012 

.4773 

.3784 


5958 

4754 


*  In  dejprees 


•x-K-  In  Locatiofi 


1,  0  =  90" 


SHEAR  FRONT  <t>,= '56.4° 


160” 


150®  140® 

INELASTIC  front  ^"120.1 


FIG. 190  RESULT  IN  RANGE  10 


results  indicate  that  the  major  stress  and  velocity  changes  occur  at  fronts,  an 
initial  front  where  the  normal  stress  rises,  which  may  be  followed  by  a  shear 
front.  The  continuous  inelastic  regions  which  occur  are  quite  narrow,  only  a  few 
degrees,  and  produce  chemges  in  the  principal  stresses  and  their  direction  which 
are  usually  quite  email  compared  to  those  at  the  fronts.,  It  is  extremely  important 
that  a  similar  behavior  can  be  expected  if  the  surface  pressure  is  not  a  step,  but 
decays  as  shown  in  Fig.  1.  This  expectation  makes  it  possible  to  solve  the  problem 
of  a  decaying  surface  pressure  approximately  by  disregarding  the  inelastic  regions 
as  being  of  secondary  importance,  but  allowing  elastic  and  inelastic  fronts  of 
discontinuity  in  the  appropriate  locations.  Because  of  the  con5>lexity  of  solving 
the  problem  with  decaying  pressure  exactly,  it  is  intended  to  utilize  the  above 
approximate  formulation  in  future  work. 
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APPENDIX  A 


Proof  for  the  Requlreiaent  T  ■»  0  at  a  Plane  Discontinuity 


In  Section  Ilb  it  was  shown  that,  for  certain  values  of  or  and  v,  a  dis¬ 
continuous  inelastic  plane  pressure  front  may  exist.  This  was  shown  on  the 
premise  that  the  principal  stress  is  normal  to  the  front,  which  is  equivalent 
to  stating  that  the  shear  y,  parallel  to  the  front  vanishes.  It  will  be  proved 
here  that  no  plane  discontinuity  with  inelastic  deformation  can  propagate  unless 
T  «  0. 


Let  y  again  be  the  direction  of  propagation  and  Xqs*  (lU)-(19)  become 


3  W  ^  Bt 
Bt 

P  ?t 


(150) 

(151) 

(152) 
(153) 
(15^) 


where  X  and  K  are  defined  by  Bqs.  (59)#  (6o).  Noting  s  »  s 

4 

yield  condition,  Eq.  (4),  becomes 

T®  +  I  >y  -  Ot^jf  -  0 


I 


(155) 


Steady-state  solutions  require 


V(C) 


.  agf(C) 


(156) 

(157) 
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V  «  a^fCC) 


(158) 


T  -  aj^f(C) 

(159) 

u  «  ^5^(0) 

(160) 

vhere  C  ■  y  “  ct. 

Substitution  of  Elqs.  (159);  (I60)  into  Eqs.  (153);  (15^)  yields,  after 
simplification 


.  5 


(161) 

1 

1 

(162) 

f  being  subject  to  the  inequalities  (64),  (65),  neither  f  nor  f  can  be  zero 
and  Eqs.  (I61),  (I62)  permit  nonvanishing  solutions  for  aj^  and  a^  only  if 


c 


2 


(163) 


X  -  0 


(164) 


The  requirement  X  =  0  violates  the  basic  condition  X  >  0  at  locations  of 
inelastic  deformation.  Therefore,  a^  and  a^  vanish,  i.e.  no  discontinuity  in 
shear  can  occur. 


Having  demonstrated  that  an  inelastic  discontinuity  in  shear  is  impossible, 
it  remains  to  be  shovm  that  the  discontinuity  in  the  normal  stress  can  not  occur 
in  a  region  with  shear,  “  /  0,  even  if  t  is  continuous.  In  this  case  Eqs.  (I50)- 
(152),  (155)  aJid  (156)“  (158)  apply  while 

T  ■  g(y;t)  (165) 

u  ■  g(y;t)  (166) 

where  g(y;t)  /  f(y  -  ct).  Substituting  Eqs.  (I65)  and  (I56),  (157)  into 
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Eq.  (155)  gives 

g®(y,t)  +  (|  -  or^i^)  f®(y  -  Ot)  -  0  (167) 

f  and  g  being  different,  nonvanishing  functions,  this  equation  can  not  be 
satisfied. 
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APPENDIX  B  -  Steady-State  Solution  in  an  Elastic  Half -Space 


As  background  for  Section  III,  the  details  of  the  solution  of  the  steady- 
state  problem  for  an  elastic  half -space  are  derived.  The  calues  of  the  stresses 
in  Cartesian  coordinates  could  be  obtained  by  integration  from  Ref.  fl]  and  the 
desired  principal  stresses  could  be  computed.  However,  it  is  just  as  easy  to 
obtain  the  latter  directly  from  the  knowledge  of  the  location  of  the  shock  fronts 
epp  and  cpg  in  Fig.  A-1,  coupled  with  the  necessity  of  uniform  stresses  for 
'Pg  >  9  >  9p  and  n  >  <p  >  9g  .  The  values  epp  and  tpg  depend  on  the  velocities  of 
the  fronts  and  ore  given  by  Eqs.  (Il4)  and  (115)* 

Designating  the  principal  stresses  in  the  region  (pg  >  9  >  epp  by  and 

Co  «  <7  it  follows  from  Eq.  (91)  that 
j  z 


c^  =  Ac 


(168) 


where  the  jump  Ac  remains  to  be  determined.  The  direction  of  c^  makes  an  angle 
(9g  -  9p)  with  the  normal  N  to  the  S-front.  The  normal  stress  c^j  ,  and  the 
tangential  stress  c^  with  respect  to  the  S-front  in  the  x-y  plane  can  be  expressed 
by  the  principal  stresses  Cj^  and  Cg  , 


Cjj  =  Ac  [co82(9g  -  <Pp)  +  8in2(9g  -  9p)] 

=  Ac  [8ln^(<Ps  “  'Pp)  +  co8^(9g  -  9p)] 


(169) 


(170) 


In  the  region  tt  >  9  >  9g  the  principal  stresses  are  c^  ,  Cg  and  c^  s  . 
The  surface  condition  requires  that  C^  =  -  1  be  vertical,  making  an  angle 
(tt  “  9g)  with  the  normal  to  the  shear  front.  The  normal  and  tangential  stresses 
(with  respect  to  the  S-front)  are  therefore 


AAW4  —  W  • 

3  z 
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FIG  A~l 


F  I  G  .  A-2 


F I G  A-3 
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where 


2  2 

cos  tPg  +  R  sin  cpg 

2  2  1 

sin  9g  +  R  cos  (pg 


(171) 

(172) 


(173) 


There  being  no  discontinuity  in  the  normal  and  tangential  stresses  at  a 

shear  front,  Ojj  and  in  Eqs.  (l69)-(l72)  can  be  equated  and  give  two  simultaneous 

equations  for  Ao  and  R.  The  stresses  and  o^  in  the  z  direction  must  also  be 

equal,  CT^  ■  a-  •  Using  the  abbreviation 
•i  •j 

2  2 

N  =  cos  9g  +  (l-2v)  cos  (<pg  “  cpp)  -  1+v  (17^) 

the  discontinuity  at  the  P-front  for  a  unit  surface  load  becomes 

Aa  -  -  cos  2tpg  (175) 


In  the  region  cp  >  the  principal  stresses  are 


=  -  1 

(176) 

cos  atpcj 

=  1  - 

N 

(177) 

CO 

to 

o 

o 

II 

(178) 

while  the  invariants  become 


,  ,  .  (l-wv2)  2  „ 

Jg  =  1  -  — ij—  +  cos  apg 


(179) 

(180) 


It  is  noted  that  cpp  and  cpg  are  functions  of  v  in  such  fashion  that  K  for  V/cp  >  1 
is  necessarily  positive,  so  that  the  condition  <  0  is  always  satisfied. 
However,  the  yield  inequality  gives  a  condition  on  o,  Eq.  (93)  in  the  text. 
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APPENDIX  C  *  Analysis  for  Range  Ib 


In  this  range  discontinuoua  elastic  stress  changes  occur  at  the  P-  and  S  fronts 
so  that  the  coobined  effect  satisfies  the  yield  condition  at  f 
following  derives  the  required  details  of  the  state  of  stress  for  f  > 


Using  an  approach  similar  to  that  in  AppendJ.x  B,  the  stresses  in  the  region 
9g  ^  ^  f  given  by  Eqs*  (l68),  and  the  normal  and  tangential  stresses 

with  respect  to  the  S-front  by  Eq.a«  (l69)>  (170).  The  principal  stress  for 
f  »  will,  in  this  range,  make  an  unknown  angle  i  with  the  normal  N  to  the 

S-front,  Fig.  A-2,  and  and  becaac  alternatively 


2  2 

"  ®1  ® 


2  2 


(181) 

(182) 


Equating  Eqs.  (I69)  to  (I8I),  and  (I70)  to  (182),  gives  two  equations  for  the 
four  unknowns  Ac,  ,  Cg  and  6,  while  the  yield  relation,  F  »  0,  furnishes  a 
third  equation.  The  three  equations  are  homogeneous  in  Ac  and  Cj^  ,  Cg  so  that 
6,  y,  •  the  stress  ratios  ^  ,  and  their  equivalent  A  can  be  conputed.  One  finds 


P  ■ 


-2/ 


12Qf‘^l  -  1 


(183) 


where  the  positive  root  is  to  be  used.  (The  negative  root  correspondo  only  to  a 
trivial  Interchange  between  and  fig  . ) 


The  principal  stress  deviator  and  the  invariant  are 


(184) 

•^1  * 

(185) 
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where  Aa  is  the  as  yet  arbitrai.'y  stress  discontinuity  at  cpp  ,  while  6  is 
obtained  from  the  equation 

cog  ■  p  cog  2(«pp  -  (pg)  (186) 

Excluding  the  trivial  addition  of  multiples  of  TT,  there  are  two  roots 
+  |6[  such  that  there  are  two  possible  values,  each,  for  y  9: 

Y»|+|fi|  (187) 

®  ■  fpg  “  ^  ±  hi  (188) 

The  value  |6|  in  the  Range  lb,  considered  here,  has  bounds  which  can  be 
established  by  the  following  reasoning.  While  the  stresses  and  at  the 
S-front  are  equal  for  flJid  the  shear  stresses  t,  t  for 

(  + ) 

respectively,  are  different.  The  condition  F  <  0,  P  *  0,  applying  for  (pg 
(pg”^,  respectively  require  that  Irj  >  jt),  so  that  the  angle  between  and 
the  normal  N  must  be  larger  than  the  one  between  imd  N,  or 

h!  >  (9s  -  9p)  (189) 

Further,  Range  Ib  is  by  definition  a  range  in  which  the  entirely  elastic 
solution  does  not  apply.  If  jsj  is  larger  than,  or  equal  to  rr  -  cpg  one  could 
select  a  shear  jTg|  <  {tj  such  that  6  *  ~  .  Because  of  the  inequality  on  , 
the  condition  F  <  0  would  apply  for  this  state  of  stress,  giving  an  elastic 
solution,  which  contradicts  the  definition  of  the  range.  The  angle  |6|  is 
therefore  bounded. 

^  -  9s  >  Ul  >  (9s  “  9p)  (190) 
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APPENDIX  D  -  Analysis  for  Reunge  Ila 


Solutions  in  this  range  have  an  Inelastic  pressure  front,  but  are  other- 
^fise  entirely  elastic. 

The  stresses  in  the  region  <  9  <  f 3  i  Fig*  A-3,  may  be  cbtedned  from 
Section  lib.  Let  be  the  as  yet  unknown  discontinuity  at  then  one  finds 


where 


(191) 


1  +  acr  /3 


(192) 


In  the  region  f  the  principal  stress  «  -  1  is  vertical;  the  stress 

g^  must  equal  vhile  Og  remains  to  be  determined,  or 


®i  “  "  ^  <^2  *  ^®'l  “  "  ^  ^3  * 

where  R  is  unknown.  (The  possibility  of  (Jg  »  -  1  being  vertical  would  be  a 
trivial  interchange  of  subscripts.) 


At  the  shear  front,  the  normal  and  tangential  stresses  and  g„  must  be 

ri  T 

continuous,  which  gives  two  equations  to  determine  the  unknown  quantities  Aff 
and  R, 

cos  ap 

Aff  - — -g  "-. - — — ::: — — —  (19^) 

(i-ft)  COB  («  -  «g)  +  (i+R)  cos  -  1 
R  -  -  1  -  (l+R)  Ac  (195) 

f  *f  ) 

To  check  the  condition  F  <  0  for  tt  >  <p  >  <pg  '  the  invariants  can  now  be 
determined  using  Eqs.  (193),  (195). 
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(156) 


or, 


-  (1  +  ?i)  Acr 

Jg  -  1  +  (1  +  R)  A©  +  {t,af  I 

After  maiiipulations  the  condition  P  <  0  may  be  brought  into  the  form 
sin  (^pg  -  29)  sin  29  <  0 
due  to  sin  2^  <  0, 


(197) 


sin  (Itqjg  -  29)  >  0 


(198) 
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APFEWDIX  g  •»  Analysis  for  Range  III 


It  van  concluded  in  Section  III  that  eolutions  without  initial  die- 
continuities  exist  in  Ran^  III.  Such  solutions  start  at  f  «  9  with  initial 
values  $  »  3>  y  2  >  which  the  differential  equations  become  singular 
so  that  their  solution  requires  special  treatB»nt.  To  obtain  aByn5)totic 
solutions  near  the  singulai'ity,  the  variables  |  and  y  are  replaced  by 
f,  A  and  respectively,  defined  in  Eqs.  (138)»  The  new  variables  are  deemed 
to  be  small  quantities,  so  that  approximate  equations  can  be  obtained  by 
retaining  in  each  ej^jresaion  only  the  leading  terms  in  the  above  quantities. 
However,  the  relative  magnitudes  of  the  three  quantities  are  not  known  before¬ 
hand,  requiring  the  retention  of  the  leading  terms  in  each  of  the  variables. 
The  determinantal  equation  (40)  becomes 


b/  - 1? .  bj. 


where 


bg  - 


12  {3  -  (I-2X)  [1  -  2V  -  IWf  /3  (1+v)]} 

(I-2X)  {i-2X(l-2v)  +  «  (1+v)  [4  +  or  (2-3^^)]} 


(199) 


(200) 


24  [1  -  2v  -f  6er*'(l-fv)) 


{l-2i(l-2v)  +  «f  (1+v)  (4  +  or  /3  (2-3i()]} 


1-fv 


3^i‘*2v; 


1 


(201) 


V  (202) 
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1 


* 


While  bp  is  positive  everywhere,  bg  is  positive  in  Range  III,  considered  here. 
Using  Eqs.  (44)-(li6)  expressions  for  3  and  y  can  be  formed  and,  after  changing 
to  the  new  variables,  lead  to 


dc 


*  M  Gli  4 
A^T)  +  1 


dc 


[B^A  +  Bg#] 


where 


.  .  1  f-  3  *  (l-2t)  [1  •  2v  -  to  /3  (Itv)l 

^  3(1“2X)  ^  1  +  (l-2v)  (I-2X) 

B  .  ■*  fl  -  2v  -  »  /3  (Itv)] 

^  ’  I-  1  +  (l-2v)  (l-2ii)  •* 

B  .  -  fs  ^  a-2x)  [1  -  2V  -  a»/3  (l.v)ll 

®  (1-2X)  I-  1  +  (l-2v)  (1-2X)  J 

where  X  is  given  by  Eq.  (202). 


(203) 

(204) 

(205) 

(206) 

(207) 


The  knowledge  of  the  nondimensional  stress  variable  A,  equivalent  to  P, 
is  not  sufficient  to  find  the  stresses,  and  one  additional  relation  is  required. 
The  most  suitable  one  is  obtained  by  adding  Eqs.  (44)  and  (45),  leading  to  an 
equation  for  (s^^  +  Sg), 


+  ig)]  »  C 


1  V 


(208) 


where 


0^  .  .  (l«) 


•• 

1 

u-2v/  Vl“2v) 

fTTe? 

m)] 

[1  +  (l-2v)(l-2X) 

(209) 


When  solving  the  three  equations  (199 )#  (203);  (2C4)  in  the  three  unknowns 
T),  A  and  OL/V,  the  first  two  are  small  quantities,  while  GL/V  must  go  to  infinity 
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in  the  limit  t  0.  (The  possibility  of  finite  values  for  this  limit  has  been 
previously  eliminated  in  Section  III  as  permitting  only  trivial  solutions 

Sj  «  ■  0. ) 

bo  and  bo  being  positive,  Eq.  (199)  is  hyperbolic  in  character,  and  permits 
o  9 

two  types  of  asymptotic  solutions.  In  solutions  of  Type  A,  “T]  and  A  are  proportional 
to  /e  , 


7)  -  D^/i 


A  *  Dg  /t 


(210) 


while  for  solutions  of  Type  B,  T,  is  proportional  to  /i  ,  while  A  is  small  of 
higher  order. 


T)  ■  /« 


A  «  Dg  t 


(211) 


where  N  >  ^  . 


For  solutions  of  Type  A,  Eqs.  (210),  the  leading  terms  on  the  right  side 
of  Eqs.  (203),  (204)  only  are  retained,  giving 


^  i  <n 


(212) 


M  -  B  A 


(213) 


Elimination  of  Gl/V  and  substitution  of  Eqs.  (210)  leads  to  a  requirement  on 
the  coefficients,  .  This  requirement  is  not  satisfied,  so  that  solutions 

of  Type  A  are  impossible. 


To  obtain  solutions  of  Type  B,  only  the  term  boTp  on  the  left  side  of 


Eq.  (199)  is  retained,  so  that 


(214) 
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Equation  (212)  applies  again,  giving 


GL  1 
V  *  2A^« 


(215) 


This  relation  gives  the  proper  sign  for  L  and  satisfies  the  requirement  for 
singularity  of  OL/V.  To  determine  the  quantity  A  it  is  noted  that  Eq.  (213) 
would  apply  if  N  lies  in  the  range  ^  <  N  <  1  for  the  e:q)onent,  so  that  in  this 
case  again  no  solutions  can  exist.  Alternatively,  assuming  N  >  1,  substitution 
of  T)  and  OL/V  gives  a  solution  for  A  proportional  to  t,  equivalent  to  N  ■  1, 
which  is  a  contradiction.  This  leaves  solely  the  possibility  N  «  1,  for  which 
case  Eq.  {2Qh)  indeed  gives  without  further  simplification  the  solution 


B, 


(216) 


Being  proportional  to  t;  the  quantity  A  is  small  compared  to  T],  so  that  -  as 
a  first  approximation  -  the  relations 


(217) 


A  ~  0 


(218) 


may  be  used.  Substitution  of  Eq.  (215)  into  Eq.  (208)  gives  after  integration 

(s^  +  Sg)  -  (219) 

where  is  an  open  constant  of  Integration,  while  the  exponent  is 


n 


(220) 


Equation  (2l8),  stating  A  ~  0,  xmplies  3  ~  3>  such  that  the  ratios  of  the 
stresses  must  be  the  same  as  at  the  inelastic  shock  front 


B 


2 


1 

2  ®1 


(221) 
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(222) 


J,  ■  8, 

1  2cr  1 


»  I  1  + 


YiF/  ^ 


(223) 


indicating  that  all  stresses  are  proportional  to  t*'.  It  is  important  that 
this  e3q)onent;  while  always  positive,  is  less  than  unity  and  usually  a  very 
small  nuaber,  of  the  order  of  1/100  (for  the  specific  case  v  ■  0,  cr  ■  0.05  one 
finds  n  ■  0.00598)*  Tiie  derivative  of  the  stresses  with  respect  to  the  angle  (p 
is  infinite  for  f  -♦  0,  and  the  small  value  of  n  indicates  a  very  rapid  stress 
rise  adjacent  to  the  singularity. 
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164.03 

163.43 

162.99 

0.5979 

0.4662 

0.3622 

0.6319. 

O.509S 

0.4118 

76.76 

78.76 

80.16 

0.9441 

0.9410 

0.9409 

0.5797 

0.4584 

0.3633 

65.03 

63.03 

61.63 

1/4 

163.22 

0.10 

0.15 

0.20 

152.10 

150.68 

150.08 

164.06 

164.04 

163.95 

164.97 
164 .62 

164.38 

0.6053 

0.4776 

0.3760 

0.6233 

0.4969 

0.3965 

84.34 

85.76 

86.36 

0.9791 

0.9799 

0.9801 

0.6012 

0.4773 

0.3784 

62.10 

60.68 

60, 03 

1 

0 

171.87 

0.10 

0.15 

0.20 

171.31 

170.52 

169.89 

171.88 

171.89 
171.89 

172.54 

172.30 

172.18 

0.6105 

0.4823 

0.3805 

0.6177 

0.4920 

0.3918 

82.43 

83.22 

83.85 

0.9703 

0.9761 

0.9789 

0.5958 

0.4754 

0.3779 

81.31 

80.52 

79.89 

1/8 

172.48 

0.10 

0.15 

0.20 

170.28 

169.55 

169.04 

172.51 

172.51 

172.51 

172.78 

172.69 

172.64 

0.6111 

0.4834 

0.3820 

0.6171 

0.4908 

0.3903 

84.67 

85.41 

85.91 

0.9887 

0.9889 

0.9893 

0.6071 

0.4817 

0.3820 

80.28 

79.55 

79.04 

1/4 

173.37 

0.10 

0.15 

0.20 

169.21 

168.70 

168.49 

173-42 

173.42 

173.41 

173  oO 

173.47 

173.  5 

0.6125 

0.4854 

0.3842 

O.6i4o 

0.6156 

0.4883 

0.3870 
0.61  III 

87.53 
88.03 
88.25 
8  >.91 

0.9960 

0.9961 

0.9961 

0.6116 

0.4852 

0.3846 

79.21 

78,70 

78.49 

1/3 

174.26 

0.10 

168. 61 

174.33 

174.33 

0.9999 

0.6140 

78.(1 

'*  In  de£?'ees  **  In  Uiia  iocv^tion  =  -  1,  0  =  9-3° 


